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Abstract 

In [9] and [12] we have introduced and studied a new paradigm for car- 
dinal interpolation which is related to the theory of multivariate polysplines. 
In the present paper we show that this is related to a new sampling 
paradigm in the multivariate case, whereas we obtain a Shannon type 
function S (x) and the following Shannon type formula: 
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This formula relies upon infinitely many Shannon type formulas for the 
exponential splines arising from the radial part of the polyharmonic op- 
erator A p for fixed p > 1. 
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1 Introduction 

The classical Shannon-Kotelnikov-Whittaker formula represents a function by 
means of the following series 



f(t)= l sw a(t _ jT) 



j = -co 

where T = — and / is a continuous function in L 2 (R) having a Fourier transform 
with support in [—n, n] . 

There have been many generalizations of this formula, see [7], [3], and ref- 
erences therein. There exists intimate relation between sampling theorems and 
Wavelet Analysis which have been studied exhaustively in [22], see also the 
references therein. 
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In the present paper we will consider a multivariate generalization of the 
Shannon theory which is based on a semi-orthogonal Wavelet Analysis using 
polysplines, and which has been recently developed in [S], [TU], PQ, [H], [TT]. In 
the case of what we call " spherical sampling" our approach provides a formula for 
the recovery of a function / from a "band-limited class PVq" by taking its values 
over the spheres centered at the origin of 1" and having radii e 3 for all j G Z. 
In the second case, of what we call "parallel sampling" we recover a function / 
from a "band-limited class PVq" by taking its values on all hyperplanes in R n 
defined by {x\ = j} for j G Z. 

An interesting feature of our results is that they use essentially all advances in 
the one-dimensional sampling theorems for Riesz basis as developed by Gilbert 
Walter in his book [22] , see also the more recent paper [23] . 



2 Shannon- Walter sampling with exponential splines 
(L— splines) 

In the present Section we will construct a formula of Shannon type which is 
based on exponential splines (these splines are called sometimes L— splines), cf. 



2.1 Shannon sampling for Riesz basis according to Gilbert 
Walter 

First, we will provide the construction of Gilbert Walter of sampling and the 
corresponding Shannon- Walter functions for semi-orthogonal Wavelet Analysis 
generated by a scaling function from a Riesz basis; the proofs follow the proofs 
outlined in [55] for the case of orthogonal scaling functions. 

Let us assume that the real- valued function <f> (t) is continuous and has 
shifts {(f> (t — which represent a Riesz basis for the Hilbert subspace Vq 

in L 2 (M) , i.e. 

Vb := clos L2 (M) {<t> (t ~ j) ■ j e Z} , 
and there exist two constants A, B > such that for arbitrary constants Cj 
holds 



A E 



k-i 2 < 



< 



2 

71 i 



following the tradition we call <f> scaling function of Vq. 
For every element / G Vq we have 

oo 

/(*)= E fMt-j)' 

7=-oo 

Assume that the space Vq has the property that for every / G Vq with 

/ (j) = for j G Z 



(1) 
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follows / = 0. Then there should exist a basis {i/jj (t)} jeZ of Vb such that a 
formula of Shannon type holds, i.e. for every / € Vo 

oo 

/(*)= E fovtER. 

j=—oo 

The problem of finding the proper conditions on 4> and the basis {tpj (t)} has 
been resolved by Gilbert Walter, [22]. Below we provide in detail his construc- 
tion. 

It is well-known that to the Riesz basis {4> (t — there exists unique 

dual Riesz basis |</> (t — j)j where duality means 

(ct>{x-j),${x-l))=8 jt , (2) 

cf. [2T], [55]. Let us assume that the functions <fi and 4> satisfy the following 
asymptotic conditions: 

0(*) = o(itr 1 - £ ), j(t)=o(\t\- i - s ) (3) 

as t — ► ±oo, tgE, 

or, which is equivalent, there exist positive constants G\ > and Ci > such 
that for all t G R holds 

l<M*)l < — 

(l + |t|) 1+£ 



'(*) 



Then a direct estimate shows that the function 

q( x ,y) = ^2<fi( x - j) <Mz/ - i) ( 4 ) 



p2 



is uniformly convergent on every compact set in 

On the other hand for every fixed x\ € R (or uniformly for from a compact 
subset in R ) we have the estimate 

C 2 

\q(xi,y)\ < 1^ (?/ ^ J') (y - ^)l T—— .,,i+e ' TTT+F- 

i, k&L (i + Fi-Jl) (l + |xi-fc|) 

We have also 

\<t>(y-j)4>(y-k)\dy< f — — — Cl dy 

Jw.{l + \y-j\) {l + \y-k\) 

< Ci f T^dy 



2 f 1 

= C q / 7—. — dy < 00. 
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The last implies the convergence of the integral 

\q(x 1 ,y)\ 2 dy. 
Thus for every element / 6 L^ (M) the integrals 



make sense. In a similar way follows that q{x,y\) £ L2 (K) uniformly for y\ 
from a compact subset in R. 

Thus, from ([T]) and the representation for the dual basis we obtain the fol- 
lowing reproduction property for every / £ Vq, 

q ix, y) f (y) dy = f (x) , (5) 

Q(x,v) f {x)dx = f (y) . (6) 

Let us note that by the definition it is not clear whether the function q is 
symmetric. 

Let us remark that from ((4]) it is obvious that 

q{t-j,0)=q(t,j). (7) 

We summarize the main results of Gilbert Walter's approach from [52] in 
the following Proposition. 

Proposition 1 Let <f>(t) be a scaling function satisfying (0), i.e. <f>(t) and its 
dual defined by satisfy <ft(t) ,<ft (t) = O (\t\ 1 as t — > ±00 for t € M. Let 
us define the function (f>* by 

We assume that the function <f>* (£) satisfies the non-zero condition 

0* (0 ^ for (el. (9) 

Then the following hold: 

1. The system of functions {q(t,j)}-^ z = {q (t — j, 0)}- eZ represents a 
Riesz basis ofVo, where q is defined in 01) 

2. The unique dual Riesz basis {Sj (i)}j- £Z corresponding to {q (t — j, 0)} jeZ 
satisfies 

S j (t) = S (t-j). 
Hence {So (f — j)}jez * s a R^ esz basis as well. 
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3. The following Shannon type formula holds for all f in the space Vq: 
/(*)=£/ 0") So for teR. (10) 



4- The function Sq (t) has a Fourier transform which satisfies 

HO 



So® 



r (0 



for £ e 



(11) 



Remark 2 We will call the function Sq Shannon- Walter function. Formula 
H10\) may be considered as the analog to the Shannon- Kotelnikov- Whittaker for- 
mula whereby the condition for band-limitedness of f (i.e. the compactness of 
the support of the Fourier transform f (£) ) is replaced by f £ Vq. In the context 
of Wavelet Analysis the frequency variable £ of Fourier Analysis is replaced by 
the index j . 

Proof. 1. We will prove that the linear operator T which is defined by the 
map 



T 



4>{'-j) (x) = q{x,j) 



has a continuous extension to Vq. For this purpose we will consider the space 
Vo of the Fourier transforms of the functions in Vq. 

Let us consider the operator T\ defined on the space Vq as the operator of 
multiplication by <f>* (£), i.e. 



Tx[f] (0 :=/(0^* (0- 

Since by the function <fi has a fast decay it follows that </>* (£) is a bounded, 
continuous and periodic function and satisfies 

|</>* (£)| > c > for£eK. 
Hence, the linear operator T\ is bounded and its inverse defined by multiplication 
with {^)* (£)^j is also bounded. 
Let us see that 

T = F~ x TiF 

where F denotes the Fourier transform and F^ 1 its inverse. Indeed, it is enough 
to check this on the elements <f> (x — j) . We have 



F 
TiF 



•(•-j)J (0=e" y ^(0 



0(0 



•(0 • 
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and taking F^ 1 shows finally that F~ 1 T 1 F <j) (• - j) (t) =q(t- j, 0) . 

So we have proved that T is a bounded invertible operator, and by a basic 
result in [5T] (p. 30 ) it follows that the system {q (t, j)}j eZ is a Riesz basis of 
V . 

2. By a theorem in [21] there is a unique Riesz basis Sj (t) which is biorthog- 
onal to {q(t-j,0)} j&z , i.e. 

(q(t-j,0),S k (t)) = 6 jk . 

On the other hand due to a change of variable r = t — k we have 

(q(t-j,0),S Q (t-k)) = f q(t-j,0)-S Q {t-k)dt 



q(T + k-j,0)-S (r)dr 
= h-j,o 

which shows that the family of functions {So(t — j)}j^ z is biorthogonal to 
{l (t ~ h 0)}jez ■ Hence, by the uniqueness it follows Sj (t) = So (t — j) . 

3. Equation (Tl0|) is obtained by the expansion of / (t) in the basis {So (t — j)}j e % 
Indeed, from 

f(t) = ^2fjS (t-j) 

j 

and by the biorthogonality relation ([6|) and ([7|) follows 

/0') = (/,9(*-i,o)> = / i . 

4. For proving (jTTJ) we need to apply the Shannon expansion (fTQ|) to the 
function </> and to take the Fourier transform. 



2.2 The exponential splines (L— splines) 
2.2.1 Preliminaries on exponential splines 

At the beginning of Wavelet Analysis the splines played and important role (see 
some history in [3] and [13]). In 1991 Chui and Wang have constructed the 
compactly supported spline wavelets, cf. [i2|. A characteristic feature of this 
Wavelet Analysis is that it is semi-orthogonal. Further Wavelet Analysis using 
exponential splines has been initiated in [6] and has been discussed in detail 
in [9], [10]. We will refer to the monograph [9] for all details of the outline on 
exponential splines following below. 

The exponential splines are defined by means of ordinary differential opera- 
tors with constant coefficients given by polynomials 

N 

L(z):=U(z-Xj), (12) 

3=1 



G 



where Xj are some real constants. For simplicity sake we will denote by A the 
non-ordered vector of the numbers Xj 

A:=[Ai,A 2 ,... ) Aj V ] ) (13) 

where some of the Xj's may repeat (the number of repetitions of Xj is the 
multiplicity [ij of this Xj ); this notation is a convenient way to avoid every 
time explicitly writing the multiplicities of the A/s. We will have the operator 

<l)-(!)-n(!-.) 

and the space Un of analytic solutions defined by 

U N = U N (A) := {/ (t) : L A (J^j f (t) = for t G m| (15) 

which is well known to be of dimension N, cf. [15] . The space Un is generated 
by the exponents 

t s e xt for s = 0,1,..., MA -1 
where A G A and /ia is its multiplicity, where obviously we have 



and the sum is over all different values A in A. 

The classical polynomial case is obtained when Xj — for all j = 1, 2, JV. 

The basic objects which we will consider are the so-called cardinal expo- 
nential splines {L— splines ) generated by the operator La (Jj) • Namely, we 
define the space 5a by putting 

5 A :={/:/(f)eC N - 2 (M)flC 00 (K\Z), f Uj+1) G U N } , (16) 

i.e. they are piecewise analytic solutions to the equation L\ [4r) f (t) = 0. 

We will consider the Wavelet Analysis generated by the B— splines arising 
from the operator L\ (4z) which are called sometimes TB— splines (cf. [17] - 
these are cardinal exponential splines having minimal support). As in the classi- 
cal polynomial spline theory there exists up to a factor only one such TB— spline, 
i.e. a cardinal exponential spline Qn [A] (t) — Qn (t) G iSa supported in the in- 
terval [0, N] and satisfying 

Q N (t)>0 for*G(0,A0, 

cf. 0. 0, m- 

It is convenient to define Qn (t) by means of its Fourier transform (cf. p. 
274 in [9]), namely we have 



FT (e~ x * - e~ 4 «) 
!N (0 = J ~n TT 1 fOT i e M. (17) 



n;=i «-Ai) 
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Many important quantities which we will need are defined by means of the 
TS-spline Q N {t). 

In the case of the classical polynomial splines Schoenberg introduced the 
so-called Euler-Frobenius polynomial. In order to define a generalization of the 
Euler-Frobenius polynomials introduced by Schoenberg we will consider 
the function (see Corollary 13.24, p. 235 in Micchelli [H]) 

^- l( * ;A) := 2^/ r L A W-A ^ (18) 

where the closed contour T surrounds all points A, G A but excludes all zeros 
of the function e 1_ x ■ If we put 

N N 

r (A) := J] (e> - A) , a (A) := ]J (e^ A) (19) 

then we obtain the following representation (see Corollary 13.25 p. 235 in [9]) 

U N _ 1 (x;X) = r{X)A N _ 1 (x;X) (20) 

where IlAr_i (x; A) is a polynomial of degree < N — 1 of A. The so-called Euler- 
Frobenius polynomial 

njv_i(A) :=njv-i(0;A) (21) 

has degree < TV — 2, cf. Corollary 13.25 in [9]. 

We have the following important Proposition (cf. [14] . [15] , or Theorem 
13.31 on p. 237 and Corollary 13.53 on p. 253 in [§]). 

Proposition 3 The polynomial T1n-i (A) has exactly N — 2 negative zeros. 

By means of the TB— spline wc may define the Euler-Schoenberg expo- 
nential spline, (see p. 254 in [9] ): 

oo 

<S> N - 1 (x;\) = XQn{x-3)- (22) 

j = -oo 

Note that the word " exponential" has nothing to do with the exponential meant 
in the present paper, but it is related to the following easy to check property 

3>jv-i (x + 1; A) = A$jv-i (x; A) . (23) 

We have the following important (cf. p. 255, Proposition 13.55 and Theorem 
13.56 on p. 256, Corollary 13.57 p. 256 in [9]). 



Proposition 4 1. The Euler-Schoenberg exponential spline in H2ty) satisfies the 
following relation for all x with < x < 1, 

$ N -i {x; A) = K e- ^ ^njv_i (x; A) . (24) 
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2. If the vector A is symmetric, i.e. satisfies A = —A, then 

§N-l =®N-l forallzeC, (25) 



,2 2, 

and £/ie polynomial <&n-i (0; z) 7^ satisfies 

(0; z) ^ /or all z e C with \z\ = 1. (26) 

In particular, we see that from ([23]) and (|24]) the following representation 
holds: 

f v; A> ) = Af (°; A ) ( 2? ) 



= (-l) N - 1 A-f+ 1 exp -J] A, ; n^(A). 



We will use further the following basic 

Lemma 5 Let the vector A be symmetric, i.e. satisfies A = —A. Then the 
following holds 

IHjv-i (-1)| < |njv_i (e«) I < injv-i (1)1 for £ e R. 

Proof. By Proposition [3] the polynomial II/v-i (A) has precisely N — 2 zeros 
Vj < 0, and by Proposition 01 1) and 2) they satisfy 

VjVu-i-j = 1 for j = 1, 2, JV — 2. 
Hence, we have the representation 

N-2 

U N - 1 {\)=D [] (A -«,-), 

3=1 

where 13 is a non-zero coefficient. By the reality of the zeros Vj < it follows 

N-2 

N-2 2 1 



|njv_i (e-*«) I = n |e-«-^| = |D| J] |e"«-v 
j=i i=i 



™- 2 ,2 



T-r ic ^-Uil , t r 1 — 2Vj COS £ + V? 



3=1 1 Jl i=i 

JV-2 

, , ' 1 - 2«j + v? 



|n w -i(i)| 
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In a similar way we obtain 



N— 2 

2 



|lW e ^)|>|£|n 

= 1^-1 (-1)1. 



2.3 Checking the conditions of Gilbert Walter 

Now we put for the scaling function 

4>{t)=Q N (t). (28) 

A fundamental result in [9] (Theorem 14.6 ) says that for arbitrary non-ordered 
vector A the system {<$> (t — i)}, e g represents a Riesz basis for the space Vb 
defined by means of the cardinal splines 5a as 

V Q = S A n L 2 (R) ; (29) 

the explicit values of the Riesz constants is found in [9] . 

According to (fT7|) the Fourier transform of the basic TB— spline satisfies the 
asymptotic 



A x * ( 1 



HO = Q N (0 = ]l-^- = o^ w j forC 



oo, 



and it is clear that the derivatives of <fr (£) satisfy similar asymptotic. 

As we mentioned before formula J2]), the basis {<f> (t ~ has a dual Riesz 

basis generated by a unique function cf>, i.e. the set of functions |^ (t — j)| 

generates a Riesz basis of Vo- For the dual function <f> one has the Fourier trans- 
form (see p. 356 in [§] ) 



(0 = C% ^ for i e R (30) 

$2»-i [Aj (7V;e*«) 



ith C s = exp I - X J I 



where A = [A, —A] , i.e. it is the symmetrized vector of A. 

Proposition 6 For arbitrary integer N > 1 and arbitrary non-ordered vector 
A of order N, the functions (j) and <j) defined by Q28\) and &30\) satisfy condition 
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Proof. The asymptotic of the function <j> is clear since Q n if) is a com- 
pactly supported function, in particular, for every integer m > the following 
asymptotic holds: 



4>{t) = o 



1 

w 



for t 



±oo. 



Now from Lemma [SJ applied to the vector A and to the corresponding func- 



tion $2jV-l 



A 



(N; e 4 ^) (whereby we use essentially the symmetry of the vector 



A ), and from formula ([27]) , it follows that the function $2^-1 
bounded, namely 



A 



(JV;e«) 



is 



|n 2JV _i(-i)| < 



3>2iV-l 



A 



(iV;e l «) < |n 2 jv-i (1)| for£e 



here the Euler polynomial H 2 n-i (A) corresponds to the vector A. 

On the other hand, we have the following equality proved by integration by 
parts: 



Z7T 
1 

~ 2^ 
1 



— 00 
00 



1 



1 



2tt 1 + 1 2 



Let us estimate the last integral. Since <J> 2 jv-i A (TV; A) is a polynomial in 

A, it follows that the derivatives of &2N-1 A (iV; e^) with respect to £ have nice 
asymptotic for £ — ► ±00. Now by the above remarks about the asymptotic of 
the derivatives of the function </>(£) and the representation formula (f3"0)) follows 
the inequality 



^W(l+i 2 ) 



< C 



for t G 



This justifies the asymptotic condition for <b in ([5]). 
■ 

For proving the asymptotic of <f> and </> we needed no restrictions on the 
vector A. However the non-zero condition |9]) is more restrictive. In fact, G. 
Walter has shown that the classical splines (this corresponds to the case of all 
Xj = in A ) generate a simple and natural Shannon type formula only in 
the case of odd degree cardinal splines, [52]. Similar is the situation with the 
exponential splines. For that reason we will restrict our attention to the case 



N = 2p. 



(31) 
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Proposition 7 1. We have the equality 

0* (0 == E * C?) e " i?j = ®N-i [A] (0, e*) /or £ e K; 

hence the non-zero condition {Pp is equivalent to the condition 

$jv_i[A](0,A)^0 /or |A| = 1. (32) 

U. For arbitrary symmetric non-ordered vector A o/ order N = 2p and for 
4> given by i28\) condition (0) holds, i.e. 

<f>* (0 ^0 for (el. 

Proof. By the definition of $jv-i in ((221) and by ([M]) follows 

E * w A ~ J = E ^ (-j) A " = t A i (°> A ) 

j'ez jez 

f-l^- 1 / N \ 
- — ex P ~E A J n JV-i(A). 



(Please note the difference between $jv-i [A] and $2iv-i A since they corre- 
spond to different vectors). Now Lemma [5] gives us the lower and upper bounds 
of IIjv-i (e^)| from which we obtain immediately ([9]), i.e. 

E ^ C?) e_i?J ^ for £ e R. 

■ 

Note that by means of Lemma [5] we have used essentially the symmetry of 
the vector A and the fact that N — 2p. We will apply the above Proposition [7] 
to the case of non-ordered vectors A = = [Ai, An] given by 

_ f \ j = -k lotj = l,...,p . . 

Ak — \\j = k for j=p + l,...,2p ' 1 j 

where fc is some integer. 

On the other hand there are other vectors A which are important for our 
multivariate theory, which are not symmetric but close to being symmetric, 
and for which the non-zero condition ([9]) holds; they have been considered in 
the paper |12j and called nearly symmetric there. For us the following vectors 
A = Afc = [Ai, Ajv] will be important: 

A f Xi +j = k + 2j for j = 0, ...,p- 1, , .s 

\ A 1+p+i = -n-fc + 2 + 2j for j = 0,...,p-l; ^ 

here k is some integer. 

We have the following analog to Proposition [71 
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Proposition 8 Let A = Afc be defined by \3$ - Then cj>* satisfies the non-zero 
condition (0). 

Proof. In the proof of Proposition 13 in the paper [T2], we proved that 

Ajv-i (0; -1) ^ 

for every k > 0. On the other hand by Proposition [3] the polynomial Hiv-i (A) 
can have only real zeros with A < 0. Since by ([20]) and Proposition 2] we have 

t> N ^ [A fc ] (0; A) = -A • s (A" 1 ) A N ^ (0; A) 

it follows that the function $n-i [Afc] (0; A) has no zeros for |A| = 1. 
■ 

Hence, in the above cases we may apply the results of G. Walter as presented 
in the previous section. By Proposition [IJ 5) it follows that there exists a 
Shannon- Walter exponential spline So £ Vq = 5a H £2 (K) such that its 
Fourier transform satisfies 

~~ 0*~(£) _ ®n-i [A] (0,e«)- (35) 
From the above Propositions we obtain the following: 

Theorem 9 Let N — 2p and A be symmetric non-ordered vector of order N, or 
A be a vector given by {3$ - Then for every cardinal exponential spline f G Vq 
there holds a Shannon type formula: 

00 

f(t)= S (t-j)f(j). (36) 



3 The multivariate case 

Our main purpose is to find a multivariate generalization of the Shannon-type 
formula (fTT))) . For that purpose we will use a Wavelet Analysis which is gen- 
uinely multivariate. We will work in the framework of the polyharmonic Wavelet 
Analysis introduced and studied in [9 . It uses piecewise polyharmonic functions 
(polysplines) which generalize the piecewise polynomial functions (and splines) 
in the one-dimensional spline theory. 

We assume that the integer p > 1 is fixed. 

We will consider the solutions to the polyharmonic equation 

A p u (x) = 

in annular domains, i.e. in annuli A a ^ — {x s K" : a < \x\ < b} given by two 
numbers a, b > 0. Here 

1=1 J 
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is the Laplace operator and A p is its p-th iterate. 

For every integer k > we will assume that we are given an orthonormal 
basis of the spherical harmonics of homogeneity degree k, which we will denote 
as usually by {Y k j , cf. [5D], [9]. Thus we have 

/ Y k ^(0)Y kita (B)d0 = 5t 1 -i a ; 

here S™ _1 is the unit sphere in R™, dO is the spherical area measure, and 5j is 
the Kronecker symbol taking on 1 for j — and elsewhere. We will work with 
the spherical variables x = rO where r = \x\ , £ §" _1 . It is known that for 
special functions of the form / (r) Y k ,i (9) one has (see [9], p. 152 ) 

A p (/(r)y M (ff)) = i^/(r).y w (0) 

where the operator Loa is given by 

L d 2 + n - 1 d k(n + k-2) 

On the other hand by the change v = logr we obtain another operator M klP 
satisfying 

^>(s)=^ M "(£) 

where 

p-l 

M Kv ( z ) ■= II ( z _ k ~ 2 ^ {z + n + k- 2- 2j) , 

3=0 

cf. [S], Theorem 10.34. Here it is important that unlike L p ^ (^) the opera- 
tor Mk, p (4-) has constant coefficients. Thus we are within the framework of 
the exponential polynomials and we may consider exponential splines for the 
operator L — M ktP defined by the vector A = = [Ai, Ajv] given by (13~4")1 . 

The following representation for polyharmonic functions in annulus is a fun- 
damental result (cf. Theorem 10.39 in [9]): 

Proposition 10 Every function u £ C°° (A a ^) satisfying A p u (x) — in the 
annulus A a j, permits the representation 

u(x) = J2 «M M Y W (8) > for x = r9 (37) 
where u k ^ (r) is a solution to 

L\ k) (■£-) u k j (r) = for a <r < b. 
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Let us assume that the integer k > is fixed. 

Further we consider the space of exponential splines generated by the op- 
erator L = Mk yP (-pi . We define the space of cardinal exponential splines and 
the Wavelet Analysis corresponding to the operator L = Mk, p (■£:) as in ([29]) , 
namely we put 

V k = SjGL 2 (R):M k , p (J^J(v) Uj+1) =0, and / e C^ 2 (R)} . (38) 

Let us recall that in the case of odd dimension n (since there are no repeating 
frequencies in A of ([34]) ) the analytic solutions of Mk, p / = are linear 
combinations of the exponentials 

j e (fc+2j)^ e (-«-fc+2+2j> j^ 1 _ 

Respectively, by the change of the variable v = log r we obtain the space 

V k :={/:/ (r) =/ (log r) , for/e^}- (39) 

Finally, we may apply Theorem^ and we will obtain the Shannon- Walter 
exponential spline given by ([33)1 which we denote further by S 1 ^' (v) . By (|36p 
Sq (v) satisfies a Shannon-type formula, i.e. for every element / G Vq^ we 
have the equality 

oo 

f(v)= E fU)Sa k) (v-j). (40) 

j=— oo 

By means of the change r — e v we have / (e v ) := / (v) and we put 

S™ (r) (logr). 
Hence, for every / £ Vq we have the Shannon-type formula 

oo 

f(r)= E f^)S^{re~i). (41) 

i=— oo 

3.1 Asymptotic of Sq for A; — ► oo 

It is essential to see what is the asymptotic of (v) and of (£) as k — ► oo. 
This will be important for the existence and the regularity of the multivariate 
Shannon type function. 
We assume as above that 

0(t) = Q N [A](t) (42) 

for the nonordered vector A = A^ = [Ai, Ajv] given by ([34]) (or ([33]) below), 
and that the non-zero condition (El) holds. 
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By formula (|55|) we have 



s { o k) (0 



(0 



$jv-i[A fe ] (0;e*«)' 
On the other hand by formula ([24]) we have the representation 



(43) 



$jv-i[A fc ] (0;A) 



(-1) 



N-l 



(44) 



(k) 

Lemma 11 The Fourier transform of the Shannon- Walter function Sq and 
its Fourier transform satisfy the following asymptotic: 
1 . uniformly for ^el 



2. uniformly for all t 6 R 

S { k) {t) = 0(1) 



for k 



for k 



Proof. By (|42|) it follows directly that for appropriate constants ko > and 
C > and for all £ £ R and k > ko we have 



(0 



nf=i w-a,) 



< c 



(£ 2 + fc 2 r 



Now we will take care of the denominator in (|43j) . We will apply a subtle 
estimate of Theorem 11 in the paper [12] (please note the conventions in the 
paper, where the size of the vector A is numbered through N + 1 and here we 
use N ) . It says that there exist two constants C\ , Ci > and an integer ko > 1 
such that for k > ko holds 

\<|^v-i(0;A)|<^ for |A| = 1. 



From relation ([20]) we obtain 

njv_!(A) = r(A)^ JV _ 1 (0;A) 

and from (|4"4"|) we obtain 

[A fc ] (0; A) = -A • s (A- 1 ) (0; A) . 

From p9[) it follows that the asymptotic of the polynomial s (A) for fc 
with A = e'* is 



S (A) 
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Hence, for k — > oo and for appropriate constant C3 > we have 

^0 (?) ~<^3 — ^ — - r- 

FFT 

2. By definition we have 

/oo 
-00 

hence by a direct estimate and a change of variable £' = £/fc we obtain 

/OO 1 

= 0;. 

■ 

Obviously the results of Lemma [TT1 hold for the Shannon function Sq^ (r) := 
Sq^ (log r) in the variable r. 

3.2 The Shannon polyspline 

We have seen that for the vector defined by (j34|) the non-zero condition (132| 

(k) 

is satisfied and we may define the Shannon- Walter exponential splines Sq . 

We will call the following function (actually it is understood in a distribu- 
tional sense) Shannon polyspline kernel: 

00 ak 

?(*0 



So (r, 9, ^):=Y,J2 S o k) to *M (9) Y kl (i/>) (45) 

00 

= J2 s o k) (r)Zk(e,^) 

00 

= J2^o h) (logr)Z k (6,ip), 

k=0 

where Z k (9,ij}) is the zonal harmonic of degree k, cf. [20]. Let us note that for 
every fixed ip the function So (r, 9, tp) is a polyspline of r9, and for every fixed 9 
the function So (r, 9, ip) is a polyspline of rip. 

We will characterize the smoothness-convergence properties of this series in 
terms of Sobolev spaces. Let us recall that the function h defined by the series 

h (x) := hk,tYk,t (9) 

k,£ 

belongs to the Sobolev space on the sphere H s (S" _1 ) for a real number s if 
and only if 

J2( 1 + k2 Y h l* < °°; 

k.t 
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cf. section 7 and section 22 in [15] . 

Lemma 12 As a function of 9 (or ip ) the Shannon polyspline Sq satisfies 
S (r, ., V) G iT (S"- 1 ) , 5 (r, 0, •) G H s (S"" 1 ) 

where 

3 

s = — n + - — e 

/or every number e > 0. 

Proof. We know that for some constant C > and for all k > 0, and 
£ = 1,2, aft holds 

see e.g. [IB], [5]. Also, and one has the estimate 
cf. e.g. [20]. 

Since the coefficients of the series ([45]) are hk,e — Sq (r) Yk,t (tp) we have 
the following inequalities 

oo 

]T (1 + fc 2 ) s hl t <J2( 1 + k 2 ) S C 2 k n - 2 < C 2 (1 + fc 2 ) S A; 2 "- 4 

k,f k,£ k=0 

which proves the theorem. 
■ 

3.3 The polyharmonic Wavelet Analysis 

In [3] we have introduced a polyharmonic Wavelet Analysis by considering 
the following basic space 

PV Q := dos L2{Rn) {/ (x) : A p f (x) = on e j < r < e>' +1 , x = r9, f G C 2p - 2 (K™)} . 

Let us note that the radii e J appear in a natural way since by the represen- 
tation ([37]) and Proposition [101 every function u G PVo has coefficients Uk,£ (r) 
satisfying 

Ufc,£ G Vtf 5 , 

where the spaces V fc are those defined in (13!)]) . Hence, we have the following 
decomposition in the sense of L 2 metric (cf. [9]) 

oo a k 
k=0 1=1 

which is generated by the expansion ([3"T]) . 
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As we have seen above, the Shannon polysplinc kernel P5|) is a function in 
some Sobolev space H s and might be a distribution in the case of a negative 
exponent. Hence, the formulas which we will write will be understood in a 
distributional sense. For two functions (distributions) / S H 8 (S™^ 1 ) and g £ 
H~ s (S n_1 ) with expansions 



jn — l 



we have the scalar product on §" 1 defined by an integral understood in a 
distributional sense 



</,<?> := / f(8)g(6)d9 = J2h. 



igk,i- 



Note that the infinite sum which gives sense to the above integral is absolutely 
convergent for / G H s (S"" 1 ) and g £ H- s (S"" 1 ) . 

In our main Theorem below we will use the spaces PVq to mimic the one- 
dimensional cardinal spline space Vq, and to formulate an analog to the Shan- 
non type formula. 

Theorem 13 For every function f S PVq Pi C(M") the following Shannon 
type formula holds: 



oo ~ 

/W>) = J s Q (^-,e^)f{eJ 



d6, 



J=-«W-i 



where So (r, 6, ip) is the Shannon polyspline function defined in and 
the integral is understood in distributional sense. 

Proof. Since by definition / £ C (R n ) it follows that for all j 6 Z we have 
/ {e^Q) € L>2, hence the following Li— convergent expansion for every fixed r > 
holds true 

oo a k oo a k 

f ^ = E E /m ( r ) Y ^ (fi) = EE /m ( lo e r ) Y ^ w - 

fc=0 fcl fc=0 fc=l 

where x = rd. 

We obtain in distributional sense the following equalities, for v = logr, 

oo p oo / oo afc \ 

£ / Sb(£,M)/(e*fl)d»= E EE^^'OM^wM 

oo ak / oo 

= EE E fM/wH^w 

/c=o£=i yi=-oo 

oo a fc / oo 

= EE E 

fc=0£=l \j = -oo 
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On the other hand since for every k > we have fki{e v ) £ Vq, by the 
Shannon formula in (1401) we obtain 



EE E ^ (v j) f M (e') F M 0/0= EE ^ ^ F ^ W 

fc=0 f=l W=-oo / fe=0 £=1 

= /(e» 
= / NO 

This ends the proof. 
■ 

Let us remark that the sampling 

{/ (k) : k £ Z n } 

is the one used in the usual multivariate generalizations of the Shannon type 
formula. In Theorem 1131 we use though the sampling 

{fk,e(e j ) :fc>0, e=l,2,...,a k , j £ Z} 

which defines a new multivariate sampling paradigm. 

4 The polysplines on parallel strips 

We consider very briefly for every integer p > 1 the Wavelet Analysis generated 
by the polysplines on strips (with periodic conditions), and the corresponding 
Shannon polyspline function and Shannon formula. 
We will introduce some convenient notations: 

x = (t, y) 6 M™ 
t£R andyeT"- 1 , 
T"- 1 := [0,27r]" _1 . 

The polyspline Wavelet Analysis on strips is defined by means of the following 
space, cf. [9], 

PV :=clos L2inn) {f(t,y):APf(t,y) = for t £ Z, / g C 2p ~ 2 (M. x T" -1 )} ; 

these are functions which are piecewise poly harmonic on every strip (J,j + 1) x 
T™ -1 , for all j £ Z , and which are 27r— periodic in the variables y. The functions 
/ g PVq permit the representation 



f(t,y)= /«(*K <K ' V> . 



k£2 
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where for every k£2" 1 the function f K (t) is a solution to the equation 

w 2 = E4 

3 

In the present case we have the vector A = Ak = [Ai, Xn] given by (|33j) 
for fc = For this A = A/- and for every integer p > 1 we obtain the operator 

Now as in (|29p we have similar definition of the cardinal exponential spline 
spaces, namely we put 

V k := j/ E L 2 (R) : L\ k) (|) f (t), yj+1) = 0, and / G C 2 ^ 2 (K)| . 

For the vector Aj. defined in by Proposition [7] condition P2"|) always 

holds. Hence, for every fc > there always exists the Shannon- Walter expo- 
nential spline Sq defined in (|35p . 

Remark 14 T/ie case /c = corresponds to the classical cardinal splines of odd 
degree considered by G. Walter in t22jj . 

Since all basic lemmata have been proved above also for this type of vectors, 
we obtain the following Shannon type formula: 

Theorem 15 For every f £ PVq n C (R n ) holds the Shannon type formula 

oo „ 

f^y')= E / S (t-j,y,y')f(j,y)dy, 

where So (t,y,y') is the Shannon polyspline distribution defined by 
S (t,y,y')= E S { K) (t)e-^e^y'). 

Again we may compare the usual multivariate sampling {/ (k) : k £ Z"} 
with our sampling paradigm 

{f K (j):KEZr-\j€Z}. 

Finally, an idea for future research might be the error analysis of the Shannon 
type formulas of Theorem Q2] and Theorem [T5j 
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